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We investigate the presence of squeezing in the weakly repulsive uniform Bose gas, in both the 
condensate mode and in the nonzero opposite-momenta mode pairs, using two different variational 
formulations. We explore the 1/(1) symmetry breaking and Goldstone's theorem in the context of 
a squeezed coherent variational wavefunction, and present the associated Ward identity. We show 
that squeezing of the condensate mode is absent at the mean field Hartree-Fock-Bogoliubov level 
and emerges as a result of fluctuations about mean field as a finite volume effect, which vanishes in 
the thermodynamic limit. On the other hand, the squeezing of the excitations about the condensate 
survives the thermodynamic limit and is interpreted in terms of density-phase variables using a 
number-conserving formulation of the interacting Bose gas. 
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I. INTRODUCTION 

The quantum optics concepts of minimum-uncertainty 
states such as coherent and squeezed states have been ap- 
plied to quantum condensed-matter systems in a variety 
of settings. The study of the Bose gas with weak repulsive 
interactions has benefited greatly from borrowing and ex- 
tending the ideas of coherent states originally developed 
by Glauber [1, 2] in optics, particularly in the context 
of the Bose-Einstein condensate [3-8] . In the traditional 
treatments of the interacting Bose gas, pairs of opposite 
momentum excitations can be created (destroyed) out 
of (into) the condensate, and one can interpret this ef- 
fect in terms of squeezing of pairs of opposite-momentum 
excitations induced by inter-particle interactions. Fur- 
thermore, squeezing within the condensate mode itself is 
another intriguing question, related to higher moments 
of the condensate-mode annihilation and creation oper- 
ators [8-15]. Thus, naively the interacting Bose gas dis- 
plays two types of squeezing super-imposed on properties 
of a coherent state that would, strictly speaking, only 
represent the behavior of a non- interacting gas. While 
"squeezing" effects have been known in some form or 
other for a long time, it is only recently that the language 
of minimum-uncertainty states has been used in descrip- 
tions of the Bose gas [9-12, 16]. Despite several studies, 
there remain significant questions concerning the exis- 
tence and physical interpretation of squeezing in the Bose 
gas ground state. In this Article, we seek insight into the 
intuitive physical meaning of squeezing in the context of 
the weakly interacting Bose gas, treating both kinds of 
squeezing mentioned above: the single-mode squeezing 
within the zero-momentum condensate mode, as well as 
the pair-wise squeezing of finite-momenta bosons. 

A coherent state encodes the physics of having a defi- 
nite phase at the expense of strict number conservation, 
which to the condensed matter community is the essence 
of Bose condensation. (See however Refs. [17-22] for 
attempts to circumvent number conservation violation.) 



As a result coherent states were appreciated early in the 
study of the Bose condensate [3-6, 8]. In addition, the 
need to incorporate ±k correlations also led to squeezing 

operators similar to exp 7cj i cL k — 7*CkC_k , which to- 
day would be called "squeeze" operators, being used for 
the interacting Bose gas since the 1960s [3, 4, 8, 13, 23]. 
In addition, some early authors have also incorporated 
single-mode squeezing explicitly in the condensate mode 
itself [8, 13-15]. 

During the resurgence of interest in the interacting 
Bose gas in the 1990s, several studies on squeezing in 
the Bose gas have been performed explicitly using the 
quantum-optics language now available. Studies of the 
quantum state of trapped condensates [11, 12] have indi- 
cated the presence of squeezing in the condensate mode 
itself, which corresponds to k = mode squeezing in 
the uniform case. Ref. [16] uses a "generalized coherent 
state" (including both single-mode and two-mode squeez- 
ing), to derive time-dependent Hartree-Fock-Bogoliubov 
equations for a non-uniform Bose gas. Refs. [9] and [10] 
have both used a wavefunction containing a squeezed co- 
herent state for the condensate mode, and the usual ±k 
pair squeezed vacua for the k^O modes, similar to our 
first variational wavefunction |sql) in Sec. II. Ref. [9] 
considers squeezing in the condensate mode (as we do 
in Sees. II and III) and focuses on regulating anomalous 
fluctuations, while Ref. [10] uses the formalism to calcu- 
late coherence functions [1, 24]. 

Our main results are as follows. We find that the con- 
densate mode is indeed squeezed, but the scaling of the 
squeeze parameter with system size is such that squeez- 
ing has no thermodynamic effects. For finite-size sys- 
tems, the presence of appreciable squeezing is determined 
by the competition of two small parameters. For Bose- 
Einstein condensates in traps, this is the same competi- 
tion that determines whether the density profile is gaus- 
sian or is given by the Thomas- Fermi approximation. We 
have also formulated the Hugenholz-Pines (H-P) theorem 
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[25, 26] in the context of our variational formulation. 
The H-P theorem enforces the absence of a gap in the 
excitation spectrum of the condensed Bose gas. (For a 
modern description, see, e.g., Ref. [27]). We use the H-P 
theorem, or the equivalent requirement of gaplessness, to 
prove that any condensate-mode squeezing present in the 
system must come from beyond a mean-field treatment 
of the theory. 

In addition, we give a physical interpretation to the 
pair-wise squeezing of boson pairs induced by conden- 
sate depletion, by using an alternate variational state, 
based on the number-conserving formulation of the Bose- 
condensed state in Ref. [21]. The finite- momenta squeez- 
ing can be expressed in a "quadrature" space of density- 
oscillation and phase operators. The squeezing of fluctu- 
ations is found to be in the density-oscillation direction. 

The paper is organized as follows. In Sec. II, we briefly 
review relevant concepts of coherent and squeezed states 
(Sec. II A), and then construct our first variational wave- 
function for the zero-temperature Bose gas (Sec. II B). 
This is used to derive the scaling of condensate-mode 
squeezing properties with system size (Sec. II C). In 
Sec. Ill we explore the manifestation of U(l) symmetry 
breaking within this formalism, formulate the relevant 
Ward identity (Hugcnholz-Pines theorem [25-27]), and 
construct the excitation spectrum. These results lead 
to additional physical inferences about the condensate- 
mode squeezing, which are presented in Sec. IIIC. In 
Sec. IV we present a second variational state, using 
density-oscillation and phase variables introduced in 
Ref. [21], and use this construction to provide a phys- 
ical interpretation of k^O squeezing. 



II. VARIATIONAL TREATMENT OF BOSE 
GAS USING SQUEEZED COHERENT 
WAVEFUNCTION 

The three-dimensional uniform Bose gas is described 
by the Hamiltonian: 



H = ^ k _ V) C|L 5 k + 4. + k 5 q-kCpCq , (1) 



2V 

k p,q,k 

where = fc 2 /2m is the free-gas dispersion, fh is the 
boson mass, and c, & are bosonic operators. The inter- 
action U is taken to be momentum-independent because 
at low enough temperatures only s-wave scattering is im- 
portant: U = 4irah 2 /rh modulo an ultraviolet renormal- 
ization term, where a is the s-wave scattering length. A 
dimensionless measure of the interaction is an 1 / 3 , where 
n = N/V is the density. 

In this Section, after a lightning review of the relevant 
coherent and squeezed state concepts (Sec. II A), we will 
introduce our first variational wavefunction |sql) and 
determine the variational parameters by minimization 
(Sec. II B), and discuss variances and squeezing prop- 
erties (Sec. II C). 




FIG. 1: Fluctuation contours of squeezed coherent states (A 
and B), squeezed vacua (C and D), and a coherent state with 
real a (E). States A and B have their squeeze phases locked 
to twice the coherence phase {cj> = 26), and are therefore 
respectively amplitude-squeezed (s < 0) and phase-squeezed 
(s > 0) states. State C has <f) — 0, s < 0. State D has larger 
s, and <f> = 7r/2, so that the tilt is n/4. 



A. Quantum States of Bosonic Systems 

Details on minimum-uncertainty quantum states can 
be found in quantum optics texts and reviews, e.g., in 
Refs. [24, 28]; we give here only a brief introduction to 
squeezed and coherent states and point out some applica- 
tions of squeezing concepts in condensed-matter systems. 
For a bosonic mode described by operators c, $ , one de- 
fines hermitian "quadrature" operators X — ^(c+c' ) and 
Y = ^(c — cr), conjugate to each other, [X, Y] = i/2, so 
that the uncertainty relation {5 2 X) 1 / 2 {5 2 Y) 1 I 2 > 1/4 is 
satisfied. Coherent and squeezed states both have mini- 
mum uncertainty. 

Coherent states have equal uncertainties in the 
quadrature directions. A coherent state can be con- 
structed by applying the displacement operator D(a) = 
exp [ac^ — a*c] on vacuum. The vacuum itself is a spe- 
cial case with a = 0. Coherent states have circular vari- 
ance profiles, centered at ((X), (Y)) = (lZe[a],Im[a]). 

Single-mode squeezed states are produced by the 
squeeze operator S = exp [76^ — 7*cc] , with 7 = 3 e 1 ^, 
whose effect is to squeeze variance profiles in the direc- 
tion indicated by <f>/2 on the quadrature plane. When 
applied to coherent states, S creates squeezed coherent 
states: 

|sq.coh) = S |coh) = SD [vac) . 

The inverted order of operators, DS [vac), is common in 
the quantum optics literature. The fluctuations of X, Y 
are the same in this alternate form, but the expectation 
values differ by the factor F — cosh(2s) + sinh(2s). The 
uncertainty contour in the X-Y plane is elliptical rather 
than circular, centered at a displaced position (a or aF). 
The uncertainties are ^e ±s along the major and minor 
axis directions. 

Applying the squeeze operator S on a vac- 
uum state produces a single-mode squeezed vacuum, 
|sq_vac(single) ) = 5|vac), which is another minimum- 
uncertainty state with distorted variances in the X-Y 
plane, but the quadrature expectations are now zero 
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(states C, D in Fig. 1). We will be more interested in 
mixed-mode squeezed vacua, 

|sq_vac(mixed)) = cxp [jc^cl - -f*c n c m ] |vac) , 

which will be the fate of k ^ states in our varia- 
tional treatment of the non-ideal Bose gas. Squeezing 
of uncertainty is now seen not in the space of the indi- 
vidual mode quadratures, (X m , Y m ) or (X n , Y n ), but 
in the mixed quadrature variables X = -^{X m + X n ), 
Y = ^(Y m +Y n ). 

After squeezed states became popular in quantum op- 
tics in the 1980s, the concept was utilized in the anal- 
ysis of several condensed matter systems. Squeezed co- 
herent states have been used to treat variationally the 
"spin-boson" model that arises in connection with dissi- 
pative tunneling [29], defect tunneling in solids, the po- 
laron problem, etc. [30-33]. Squeezed states have also 
been used for polaritons [34], exciton-phonon systems 
[35], many-body gluon states [36], bilayer quantum Hall 
systems [37], phonon systems [38], and attractive Bose 
systems on a lattice [39]. 



B. Variational wave function and minimization 

For a uniform condensate, the macroscopic occupation 
is in the zero-momentum state, so we will use a coher- 
ent occupation of the k = mode only: D = D = 

exp ac\ — a* Co , with coherence parameter a — f e 10 . 

Intuitively, a corresponds to the order parameter for Bose 
condensation. 

We will apply a mixed-mode squeeze operator for 
each opposite-momenta mode pair. Thus the variational 
ground state is [sql.gr) = 5|coh) = SD [vac), with 



^ = n ^ = n exp \ (^c-k - 7kCkc- k ) 



7k + 7-k = 2s k e^ k . 

Note that this automatically includes single-mode 
squeezing for the condensate (k = 0) mode, with squeeze 
parameter 70. Our variational wavefunction for the uni- 
form interacting condensate is thus a squeezed coherent 
state for the k = mode and a mixed-mode squeezed 
vacuum for each k^O mode pair. 

Minimization of the wavefunction locks the squeeze- 
parameter phases of each momentum-pair mode to twice 
the phase of the k = coherence parameter, i.e., 0k = 29 
for all k. In the following, we simply start with 0k = 20 
to avoid typing (0k — 29) arguments. 

To determine the variational parameters, we need 
to minimize the expectation value of the Hamiltonian 
(1). Expectation values in the variational ground state 
[sql.gr) are calculated using the relations D^c^D 



c k + aSk.o an d S'cuS — cosh(sk)c k + sinh(sk)e 



i0k At 



The required quantities are (N^) — (c^Ck) and (Hint)- 

(iVk) = sinh 2 ( Sk ) + A^ k ,o, (2) 

and 

+ ^y Nc I J2 cosh( Sq ) sinh(s q ) + 2 ^ sinh 2 (s q ) I 
I q q ) 




cosh(s q ) sinh(s q ) 



+ 2 



^sinh 2 (s q ) 



Here we have defined N c = f 2 [cosh(2s ) + sinh(2s )]- 
We can now minimize (H) = 2k( £ k — fi)(Nk) + (#mt) 
with respect to N c and Sk- This yields 



H = 2Un + Um - 2Un c 



(3) 



and 



sinh(2sk) = 



-Um 



V(e k - y. + 2Un) 2 - {Um) 2 



Here we have defined M = X)q(CqC-q) = N c + 
i sinh(2s q ) , and m = M/V. We will show in 
Sec. IIIB that the denominator appearing in sinh(2sk) 
is the quasiparticle spectrum Ek, and that \i — 2Un = 
— Um. Therefore 



sinh(2s k ) 



-Um 



-Um 



E k ^e'i + 2(Um)e k 



(4) 



C. Expectation values, variances and squeezing 



Condensate mode 
operators are (60) = 



Expectation values of the bosonic 



N c e ie and (c, 



It is 



interesting to contrast this with Bogoliubov's mean-field 
prescription, (c ) = (cj) = VNq- Since N — N c + 
sinh 2 (so), the k = squeezing parameter measures the 
deviation of our model from mean-field physics. This will 
be discussed further in Sec. IIIC. 

Since we have used a squeezed coherent state for the 
condensate mode, the expectation values and fluctuations 
of the quadrature operators Xq = i(co + cj) and Yq = 



are identical to those for a squeezed coherent 



-k- 



state in quantum optics (Sec. II A, Refs. [24, 28, 40]): 
(X ) = N c cos6, (Y ) = N c sm9, and 

(S 2 X ) = \ [e 2s °cos 2 + e- 2s °sin 2 0] , 
(5 2 %) = J [e 2so sin 2 9 + e- 2s ° cos 2 9] . 
The fluctuations along major (minor) axis directions are 

I „±s 
2 e 
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operators 



FIG. 2: Squeezing of the condensate mode. The variances are 
very much exaggerated compared to the radial distance from 
the origin. The displayed variance to radial distance ratio, if 
drawn to scale, would correspond to N ~ 0(1O 3 ), in which 
case a weakly interacting condensate (an 1 ' 3 < 0.1) would be 
barely squeezed. 



Eq. (4) shows that the squeezing parameter Sk is neg- 
ative, becomes large for small k, and diverges for k = 0. 
This infrared divergence indicates that sinh(2so) scales 
as a positive power of the system size, i.e., 

sinh(2s ) = - 7 iiV 72 and |s | ~ O(lniV) . 

The 7's can be extracted from finite-size considerations 
on Eq. (4). Noting that the lowest single-particle state in 
a box of volume V has energy eo = 3/i 2 /8mF 2 / 3 , we get 
72 = 1/3. The same result can be obtained for a power- 
law trap. The exact number 71 is probably geometry- 
dependent. However, ignoring factors like 2 and it and 
the difference between n and n C) we find 71 ~ V n 1 / 3 a; 
thus 



sinh(2s ) 



-vV/SaTV 1 / 3 . 



(5) 



upto a factor of order 1. In the true thermodynamic limit, 
V n 1 / 3 a7V 1 / 3 — * 00, so the variance profile is squeezed 
infinitesimally thin in the radial direction. The exten- 
sion of the variance profile in the phase direction, ~ 
O ([n 1 / 3 ^ 1 / 4 ^ 1 / 6 ), although diverging, is still infinitesi- 
mally small compared to the radial distance (~ yN) of 
the state from the origin in the quadrature plane. This 
is symptomatic of the fact that the squeezing has no 
thermodynamic effects, which we show more explicitly 
in Sec. IIIC. 

For finite-size systems, there is significant squeezing 
only for n'^aJV 2 / 3 3> 1. Note that this is the same con- 
dition that decides whether the Thomas-Fermi approxi- 
mation for a trapped condensate is valid or not. 

Nonzero-momentum modes — The k ^ modes in 
the wavefunction |sql) have the structure of two- mode 
squeezed vacua. The operators have zero expectation 
values, (6k) = (c k ) = 0. There is also no squeezing in 
quadrature operators defined within a single mode: the 
usual Xk, Ik have zero expectation values and equal fluc- 
tuations. 

Squeezing can be seen if one defines the mixed-mode 



X, 



k,-k 



V2 



X. 



2V2 



(ck + c_ 



at , *t 



These quadrature operators have zero expectation values 



and unequal (squeezed) variances | e ±Sk . 

The formalism of this section also allows us to calcu- 
late the occupancies iVk = (c k Ck) of the non-condensate 
modes: iV k = sinh 2 s k = 5 [(e k + Um)/Ek - 1]. This 
is consistent upto mean field order with standard treat- 
ments, e.g., Refs. [14, 27]. 



III. SYMMETRY BREAKING AND 
GOLDSTONE'S THEOREM 

We now investigate the symmetry-broken nature of the 
ground state of the Bose gas. The ground state has a par- 
ticular phase, thus spontaneously breaking a continuous 
U(l) symmetry present in the Hamiltonian. Symmetry- 
broken ground states satisfy a Ward-Takahashi identity 
reflecting the invariance of the ground-state energy un- 
der shifts of the ground state by the symmetry opera- 
tion in question. According to Goldstone's theorem, a 
phase with broken continuous symmetry should have a 
gapless mode. The Ward identity and gaplessness, both 
being consequences of the same phenomenon of spon- 
taneous symmetry breaking, are equivalent conditions 
and can generally be derived from each other. In the 
case of the Bose gas, the corresponding Ward identity is 
known as the Hugenholtz-Pines (H-P) theorem [25, 26]. 
It is the condition for gaplessness as well as a conse- 
quence of the invariance of the ground-state energy un- 
der shifts of the U(l) phase. The H-P theorem reads 
(i = S n (0,0) - £12(0,0), where £n(k,w) and £i 2 (k,w) 
are the normal and anomalous self-energies. 

In Sec. Ill A we use the fact that the Hamiltonian has 
a U(l) symmetry while the ground state (and hence our 
variational wavefunction) does not. In Sec. IIIB we cal- 
culate the excitation spectrum by constructing a single- 
quasiparticle wavefunction, and impose the requirement 
of gaplessness. The two considerations lead to the same 
condition for the variational parameters, which is com- 
forting in light of Goldstone's theorem. The condition 
should be equivalent to the H-P theorem. In Sec. IIIC 
we compare our Ward identity with the Hartree-Fock- 
Bogoliubov (mean-field) form of the H-P theorem, and 
hence evaluate the importance and effects of condensate- 
mode squeezing, sq. 



A. U(l) symmetry breaking 

In our variational wavefunction, the symmetry-broken 
nature of the ground state appears as the definite phase of 
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the coherence parameter a = \a\ e l6 . A shift of this phase 
would obviously change the wavefunction, but should not 
affect the ground-state energy, since the Hamiltonian is 
[/(l)-invariant. This requirement will give us the Ward 
identity for our formalism corresponding to the H-P the- 
orem. 

We examine the transformation 



a —* a = ae 
7k -» 7k = « k e 



fe i(e+x) ; t)^b, 



so that the ground state is shifted, | sql.gr) — > 
sql.gr) = (SY (L>y |vac). We consider infinitesimal 

A, so that a w a(l + iX), and N c = 7V C (1 + A 2 ). 
The shift in the thermodynamic potential is 



5(H) = (sql.gr H sql.gr) — ^ sql. gr| if | sql.gr 

= \ 2 N c {-n + 2Un- Un c } . 

We now use the requirement that the grand-canonical 
energy should not be changed by a shift of the ground 
state phase, i.e., 5(H) = 0. Thus we get 



fi = 2Un - Un c = 2Un - Un + U 



sinh 2 (so) 



V 



■ (6) 



Within the variational formalism, this is our equivalent 
of the H-P relation. 



B. Excitation spectrum 

We first construct a wavefunction for a Bose-gas 
ground state with a single-quasiparticle excitation added 
on top of it: 

|sql.cx(ki)) = SDc^ |vac) . 

The idea is that, since the squeeze operator S represents 
interaction effects in the present formalism, the particle 
creation operator should produce a Bogoliubov quasi- 
particle when used in conjunction with S . 

One can evaluate matrix elements of and H lnt in 
the state |sql.ex(ki) ) just as was done in the ground 
state | sql.gr). The calculation is lengthier but straight- 
forward. Calculating (H), one now obtains the excita- 
tion spectrum as the (grand-canonical) energy of the new 
state with respect to the ground state. 

£k = 

(sql.ex(ki)| H |sql.ex(ki)) — (sql.gr| H |sql.gr) 
= cosh(2si Cl ) [eki — £t + 2Un] + sinh(2sk 1 ) • Um 

= [(e kl - fj, + 2Un) 2 



For the spectrum to be gapless, one requires /j, — 2Un = 
±U m. The positive sign is inconsistent (resulting in n c = 
0). Therefore, we have 



/j, = 2Un — Um 



(7) 



Taken together with Eq. (3), this is indeed identical to 
the condition (6) obtained from consideration of U(\) 

symmetry breaking, as expected. 

The spectrum we have is thus = \J t\. + 2(Un c )e k = 
\/ e k + 2(Um)ek , which may be contrasted with the Bo- 
goliubov spectrum = y/e^ + 2(J7n )ek • 



C. Inferences on condensate mode squeezing 

In the thermodynamic limit, N — > oo, Eq. (5) im- 
plies sinh 2 (s ) w -isinh(2s ) - Vn^aN 1 / 3 . Thus, 
in our Ward identity Eq. (6), the contribution from s Q 
vanishes as 7V~ 2 / 3 for macroscopic systems. The ef- 
fect of condensate-mode squeezing on other thermody- 
namic quantities and equations similarly vanishes in the 
N — > oo, V — > oo limit, since so generally appears as 
sinh 2 (so) or sinh(2so) in equations involving extensive 
quantities. 

At the mean field Hartree-Fock-Bogoliubov (HFB) 
level, En = 2Un and S12 = Uno, so that the H- 
P theorem is ji = 2Un — Un . Comparing with our 
form jj, = 2Un — n c = 2Un — m, we conclude that the 
SD I vac) formalism reduces to the mean-field HFB re- 
sults if N = N c = M. Noting from Eq. (2) that 



No = (4 co) =N C + sinh 2 (s ) 



(8) 



(Um) 2 ] 1 ' 2 



the condition for our formalism to be restricted to mean- 
field physics is s = 0. The k = squeezing parameter 
is thus a measure of the deviation of the formalism from 
HFB physics. The point is further emphasized by rewrit- 
ing Eq. (8) as (cjc ) = (cj)(c ) + sinh 2 (s ), which shows 
that sinh 2 (so) acts as a correction to mean-field type de- 
composition. The argument can be inverted to state that, 
at mean field level, the weakly interacting T = Bose 
gas has no squeezing in the zero-momentum mode. Since 
so = at mean field level, the k — squeezing must 
come from beyond mean field. 

It may seem tempting to try to identify which diagrams 
contribute to sq = 0, i.e., to identify contributions to £12 
or £n of the form [7v / n 1 / 3 aA -2 / 3 . Note however that 
these would be non-extensive contributions, which are 
(not surprisingly) not readily found in the literature. 

Note that, in contrast to the k = mode squeezing, 
the nonzero ±k mixed-mode squeezing in the k^O modes 
is present at mean field level already. 



IV. SQUEEZING IN "FIXED-iV" EXCITATIONS 

In this Section, we will introduce and study a second 
variational formulation of the interacting Bose conden- 
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sate in order to give a more physical interpretation of the 
squeezing of the nonzero-momentum modes. The formal- 
ism will be based on the the bosonic operators introduced 
by A. E. Ruckenstein in Ref. [21]. 



A. Bosonic fields and excitation Hamiltonian 

Ref. [21] presents a current algebra approach to for- 
mulating a number-conserving description of the Bose 
condensate. The Hamiltonian is written in terms of den- 
sity and current operators h(r) = c^(r)c(r) and j = 
-^[ct(r)Vc(r)-Vct(r)c(r)]. 

The density fluctuation operator fj, defined as fj(r) = 
n(r) — (n(r)) = n(r) — n G (r), and the phase operator 
4>, defined by j = " V<ft(r), are canonically conjugate. 
Defining the linear combinations b(r), b'(r) with 



6(r) = 



1 



r)(r) + 2m G (r)0(r) 



the Hamiltonian in [21] takes the form 

H = E GS [n G (r)} + H x fn G (r), 6(r), &t(r) 



Eqs describes the mean-field ground state, and mini- 
mizing this functional gives an equivalent of the Gross- 
Pitaevskii equation which determines n G (r). In this Arti- 
cle we concentrate on the uniform case, nc(r) = n G . We 
are more interested in the excitation Hamiltonian Hx 
which describes the low-lying, large-lengthscale excita- 
tions. In momentum space, Hx reads 



H x = J^k&l&k 



k#0 

+ \Un G (Hplk + ^-k + 2 6 k S k ) , (9) 

k#0 

modulo an additive constant. Here 5 k = J r 6(r)e~* k ' r . 

The Hamiltonian (9) looks identical to that derived 
by Bogoliubov. However, the operators in the Bogoli- 
ubov Hamiltonian arc the original bosonic operators c, 
c\ rather than the peculiar bosons b, b\ that we have 
here. The interpretation is very different; the Bogoli- 
ubov picture involves an order parameter and ±k pairs 
can appear from or disappear into the condensate, while 
in the fixed-iV picture, there is no order parameter. Hx 
should not be regarded as a quasiparticle Hamiltonian, 
but rather as the Hamiltonian describing low-lying den- 
sity and current oscillations of the system at a fixed total 
particle number. It is then no surprise that Eq. (9) does 
not conserve the number of bosons, J r w(r)b(r). Our 
reason for using this formalism is that the b bosons can 
be interpreted in terms of density fluctuation and phase 
operators. 

Introducing Fourier transforms of the density fluc- 
tuation and phase operators, ?)k = J ?)(r)e~* k ' r , and 



0k = f r 4>( r )e lk r - we can express mixed-mode hermi- 
tian quadrature operators as 

*k,-k = \ (bk + S-k + bl + foL k ) = (r) k + vt) , 

*k,- k = \ { (6k + S-k -it- &l k ) = ! f(k + 4>l) ■ 

(10) 

Note that A"k,-k and Yk,_k here are different from the 
quadrature operators defined in Sec. II C because the 
bosons b, W have different meanings from c, $ . 

B. Variational treatment, squeezing 

Let us define the reference state |ref) as the vacuum 
for the 6k bosons. We now introduce the following wave- 
function as a variational state for the system: 

|sq2) = J] exp [lui>tbl k - 7k6k&-k] |ref ) = S |ref ) , 

k/0 

with the usual 7 k + 7-k = 2sk e*^ k . 

The state |ref) itself is determined from the Eqs 
part of the theory. The expectation values in our 
variational state |sq2) are (o]j>]<) = sinh 2 (sk) and 
(blb-k) = \ sinh(2s k )e-^ k . We will minimize (H x ), not 
(Hx — LiN) . This is because the number of excitations 
are not conserved, and /x is determined in the Eqs part 
of the theory. Minimization leads to 0k = 0, and 



sinh(2sk) = 



-Un G 



y/ei + 2J7n G e k 



One can also calculate the excitation spectrum from 
this alternate variational procedure. As in Sec. IIIB, 
we can construct the excited state |sq2.ex(p)) = 
S b' p | ref ). Using expectation values in states |sq2.gr) 
and |sq2.ex(p)), the dispersion relation is found to be 

E(p) = « /Cp + 2 (Uhg) e p . This is the Bogoliubov spec- 



trum, assuming n G = n Q . This demonstrates that our 
new variational formulation captures the physics of the 
weakly interacting Bose gas at least up to mean field 
level. We are therefore justified in using the formalism 
based on the state |sq2) to draw conclusions about the 
k^O mixed-mode squeezing. 

Just as in Sec. II C for the state |sql), and in Sec. II A 
for a general two-mode squeezed vacuum, the state 
|sq2) displays squeezing in the plane of mixed-mode 
quadrature operators A^-k and Yk,_k- However, at 
this stage the relevant quadrature operators are physi- 
cally meaningful: X k _ k = 4^7 (j}k + fj^j and Y k _k = 

(^ k + ^ik) ' as denned in Eq. (10). Since Sk is nega- 
tive, squeezing is along the A"k.-k direction (Fig. 3). The 
squeezing is larger for lower momentum. 
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FIG. 3: Squeezing of nonzero-momentum modes, in terms 
of "physical" variables. The variables are Fourier transforms 
of density fluctuation and phase, in units of ug and Uq 1 re- 
spectively. The density fluctuation variables are squeezed. 
The squeezing shown here is very moderate, i.e., a large- 
momentum mode or a condensate with small number of par- 
ticles. 



Thus our study of the alternate variational formula- 
tion, in terms of the bosons of the "fixed- N" theory [21], 
has allowed us to express a well-known squeezing phe- 
nomenon in terms of variables that have the very physi- 
cal meaning of density fluctuations and phases, albeit in 
momentum space. This may be regarded as a new for- 
mulation of the old idea, attributed to Feynman [41, 42], 
that in a repulsive Bose condensate, density fluctuations 
should be suppressed (or in modern language, squeezed). 



V. DISCUSSION 

In summary, we have addressed in detail the issue of 
squeezing in various modes of the ground state of a uni- 
form condensate, using two different variational wave- 
functions. 

For squeezing in the condensate mode, we have pre- 
sented a clear analysis of the scaling of squeeze parameter 
s with system size, using our first wavefunction |sql), 
resulting in the scaling relation e 2s ° <~ — V n 1 / 3 aiV 1 / 3 . 
This leads to the conclusion that while the ground state 
is indeed squeezed (with the uncertainty profile distortion 
even diverging for TV — > oo) the squeeze parameter nev- 
ertheless has no thermodynamic effects. For finite-size 
systems, such as condensates in traps, we have identified 
that the Thomas-Fermi regime (an 1 / 3 > N- 2 ^) is the 
interaction regime where one expects to see appreciable 
squeezing of the ground state. 

Our second wavefunction |sq2) is devised specifically to 
address the issue of pair squeezing in the non-condensate 
opposite-momenta mode pairs. Using results from one 
of the Z7(l)-invariant formulations of the Bose conden- 
sate [21], we have provided an interpretation of this pair 
squeezing in terms of variables representing density and 
phase excitations. 

We now make contact with relevant results in the lit- 
erature. It is worth pointing out that our treatment of 
gaplessness, where imposing the Hugenholtz-Pines theo- 
rem leads to the condition m = n c , is actually equiv- 
alent to the Popov approximation [27, 43, 44] where 



anomalous pair correlation functions (m in Ref. [43]) are 
neglected. This is a simple and direct way to imple- 
ment gaplessness. In Ref. [9], a more involved procedure 
for satisfying the Hugenholtz-Pines theorem leads to a 
macroscopic condensate-mode squeezing. However this 
contradicts the scaling relationship e 2s ° <~ — TV 1 / 3 that 
we have derived here directly from the minimization of 
variational parameters. Other scattered previous discus- 
sions of squeezing in the condensate mode [10-12, 14, 15] 
have not addressed clearly the role of this squeezing in 
the thermodynamic limit. Finally, concerning the den- 
sity fluctuation operators borrowed from Ref. [21], we 
note that similar operators have appeared in other fixed- 
TV formulations of the condensate ground state, e.g., in 
Ref. [20]. 

We end by pointing out some open problems. 

Our results on the condensate-mode squeezing prompts 
questions about the presence of squeezing in trapped con- 
densates. Study of the quantum state of trapped con- 
densates, cither experimentally through quantum state 
tomography methods or theoretically, is essential for veri- 
fying our finite-size scaling relation e 2s ° <~ — \/n 1 / 3 aiV 1 / 3 . 
Rcfs. [11, 12] have reported Q-function and Wigncr func- 
tion calculations of the condensate quantum state, show- 
ing squeezing in a number of cases. However, no sys- 
tematic analysis of the TV-dependence or interaction- 
dependence of the squeezing parameter is available. 

Another question related to the quantum state of the 
condensate mode is the possibility of non-classical fea- 
tures other than squeezing. A whole number of quan- 
tum states are studied in quantum optics (Fock, ther- 
mal, squeezed Fock, etc.) and it is intriguing to ask if, 
for example, using a squeezed Fock state instead of a 
squeezed coherent state for the k = mode would gain 
us complementary insight. Also, other quantum states 
might be helpful in exploring k ^ physics beyond the 
mean-field level physics we have extracted here for the 
non-condensate modes. Inclusion of other quantum-state 
features might also be fruitful for a variational descrip- 
tion of finite-temperature, two-dimensional, or trapped 
Bose gases. 

A real-space variational procedure using wavefunctions 
of Jastrow form has often been employed to describe in- 
teracting Bose condensates [47-49] . A natural question is 
the relation to our variational description. Presumably, 
the success of the so-called Bijl-Dingle- Jastrow wavefunc- 
tion is due to its correctly capturing correlations such as 
those we have discussed in terms of squeezing. However, 
it remains unclear how to extract from real-space Jas- 
trow wavefunctions the momentum-space squeezing pa- 
rameters of the type included in our SD state |sql). 
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